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Abstract 



We introduce twist-two quark interference fragmentation functions in he- 
licity density matrix formalism and study their physical implications. We 
show how the nucleon's transversity distribution can be probed through the 
final state interaction between two mesons (7r + 7r~, KK, or irK) produced in 
the current fragmentation region in deep inelastic scattering on a transversely 
polarized nucleon. 



Submitted to Physical Review Letters 



^Email address: jaffe@mitlns.mit.edu 
* Email address: jin@ctpa02.mit.edu 
§Email address: jtang@mitlns.mit.edu 

*This work is supported in part by funds provided by the U.S. Department of Energy (D.O.E.) 
under cooperative research agreement #DF-FC02-94ER40818 



The quark transversity distribution in the nucleon is one of the three fundamental distri- 
butions which characterize the state of quarks in the nucleon at leading twist. Measurements 
of the other two have shed considerable light upon the quark-gluon substructure of the nu- 
cleon. The transversity distribution measures the probability difference to find a quark 
polarized along versus opposite to the polarization of a nucleon polarized transversely to 
its direction of motion |l|-§]|. It is identical to the helicity difference distribution in the 
non-relativistic limit where rotations and boosts commute. However, we have learned from 
9 a/ ' 9v 7^ 5/3 and the recent measurement of the spin fraction carried by quarks in the 
nucleon, X « 0.2 |J, that the quarks inside the nucleon cannot be non-relativistic. The 
difference between the transversity and helicity distributions is a further and more detailed 
measure of the relativistic nature of the quarks inside the nucleon. 

The transversity distribution measures the correlation of quarks with opposite chiral- 
ity in the nucleon. Since hard scattering processes in QCD preserve chirality at leading 
twist, transversity is difficult to measure experimentally. For example, it is suppressed like 
0(m q /Q) in totally-inclusive deep inelastic scattering (DIS). Ways have been suggested to 
measure the transversity distribution. These include transversely polarized Drell-Yan [|TJ, 
twist-three pion production in DIS j@,|J, the so-called "Collins effect" as defined in single 
particle fragmentation J7|, and polarized A production in DIS |M. However, each of these 
has drawbacks. The Drell-Yan cross section is small and requires both transversely polar- 
ized quarks and antiquarks in the nucleon. Twist-three pion production is suppressed by 
0(1/Q), and a grfi term must be subtracted to reveal the transversity dependence. The 
Collins effect in single particle fragmentation requires a residual final state interaction phase 
between the observed particle and the rest of the jet which we believe to be unlikely (see 
below). Finally polarized A production suffers from a likely low production rate for hyperons 
in the current fragmentation region and an as yet unknown and possibly small polarization 
transfer from w-quarks to the A. 

In this Letter we develop another way to isolate the quark transversity distribution in the 
nucleon that is free from many of these shortcomings. We study semi-inclusive production 
of two mesons (e.g. 7r + 7r~, ttK, or KK) in the current fragmentation region in deep inelastic 
scattering on a transversely polarized nucleon. The possibility to measure quark transversity 
distribution in the nucleon via such a process was first suggested by Collins and collaborators 
H (see also Ref. ||). Our analysis focuses on the interference between the s- and p-wave of 
the two-meson system around the p (for pions), K* (for nK), or the <fi (for kaons). We make 
explicit use of two-meson phase shifts to characterize the interference. Such an interference 
effect allows the quark's polarization information to be carried through k + x fc_ • S±, where 
k + , k_, and S± are the three- momenta of 7r + (K), ty~ (K), and the nucleon's transverse 
spin, respectively. This effect is at the leading twist level, and the production rates for 
pions and kaons are large in the current fragmentation region. However, it would vanish by 
T-invariance in the absence of final state interactions, or by C-invariance if the two-meson 
state were an eigenstate of C-parity. Both suppressions are evaded in the p (7r + 7r~), K* 
(ttK), and 4> (KK) mass regions. 

The final state interactions of 7T7t, ttK , and KK are known in terms of meson-meson phase 
shifts. From these phase shifts we know that s- and p-wave production channels interfere 
strongly in the mass region around the p, K*, and meson resonances. Since the s- and 
p- waves have opposite C-parity, the interference provides exactly the charge conjugation 
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mixing necessary. Combining perturbative QCD, final state interaction theory, and data 
on the meson-meson phase shifts, we can relate this asymmetry to known quantities, the 
transversity distribution we seek, and to a new type of fragmentation function that describes 
the s- and p-wave interference in the process q — > 7t + 7t~(ttK, KK). Unless this fragmentation 
is anomalously small, the measurement of this asymmetry may be the most promising way 
to measure the quark transversity distribution. 

Earlier works |7| have explored angular correlations of the form ki x k 2 • S±, where ki and 
k 2 are vectors characterizing the final state in DIS. The simplest example would be k\ = q 
and ki = k n , the momentum of a pion. These asymmetries, however, require that the final 
state interaction phase between the observed hadron(s) and the rest of the hadronic final 
state must not vanish when the unobserved states are summed over. We believe this to be 
unlikely. We utilize a final state phase generated by the two-meson final state interaction, 
which is well understood theoretically and well measured experimentally. 

The results of our analysis are summarized by Eq. (|14|) where we present the asymmetry 
for 7r + 7r~ (ttK, KK) production in the current fragmentation region. Current data on nn, 
ttK, and KK phase shifts are used to estimate the magnitude of the effect as a function of 
the two- meson invariant mass (see Fig. |2|). 

We consider the semi-inclusive deep inelastic scattering process with two-pion final states 
being detected: eN± — > eV + 7r~X. The analysis to follow applies as well to 7rK or KK pro- 
duction. The nucleon target is transversely polarized with polarization vector S^. The 
electron beam is unpolarized. The four-momenta of the initial and final electron are de- 
noted by k = (E, k) and k! = (E', k'), and the nucleon's momentum is P M . The momentum 
of the virtual photon is q = k — k', and Q 2 = —q 2 = —4EE'sm 2 9/2, where 9 is the electron 
scattering angle. The electron mass will be neglected throughout. We adopt the standard 
variables in DIS, x = Q 2 /2P-q and y = P-q/P-k. The cr[(^)i=o] and p[{^^)i=l] resonances 
are produced in the current fragmentation region with momentum P^ and momentum frac- 
tion z = Ph ■ qjq 2 ■ We recognize that the tttt s-wave is not resonant in the vicinity of the p 
and our analysis does not depend on a resonance approximation. For simplicity we refer to 
the non-resonant s-wave as the "<r". In their work on the two pion system ||10|| , Collins and 
Ladinsky made the unphysical assumption of a narrow s-wave resonance interfering with a 
(real) continuum p-wave, neither of which appears in 7T7r-scattering data. 

The invariant squared mass of the two-pion system is m 2 = (k + + kJ) 2 , with k + and 
k- the momentum of 7r + and 7r~, respectively. The decay polar angle in the rest frame of 
the two-meson system is denoted by 0, and the azimuthal angle <fi is defined as the angle 
of the normal of two-pion plane with respect to the polarization vector S± of the nucleon, 
cos0 = k + x k^ ■ S±/\k + x This is the analog of the "Collins angle" defined by the 

7r + 7r~ system 0. 

To simplify our analysis we make a collinear approximation, i.e., 9 « 0, in referring the 
fragmentation coordinate system to the axis defined by the incident electron (the complete 
analysis will be published elsewhere JTTJ] ) . At SLAC, HERMES, and COMPASS energies, a 
typical value for 9 is less than 0.1. Complexities in the analysis of fragmentation turn out 
to be proportional to sin 2 ^ and can be ignored at fixed target facilities of interest. In this 
approximation the production of two pions can be viewed as a collinear process with the 
electron beam defining the common e% axis. Also we take S± along the e\ axis. 
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Since we are only interested in a result at the leading twist, we follow the helicity density 
matrix formalism developed in Refs. in which all spin dependence is summarized in 

a double helicity density matrix. We factor the process at hand into basic ingredients: the 
iV — > q distribution function, the hard partonic eq — > e'q' cross section, the q — > (a, p) 
fragmentation, and the decay (o~,p) — > ttti, all as density matrices in helicity basis: 

d 6 a 



dx dy dz dm 2 d cos 6 i 



J H'H 



sp hih[ 
J~ H'H 



d 2 o~{eq — > e'q') 


h 2 h' 2 


' d 2 M ' 


HiH[ 


' d 2 V ' 




dx dy 




dz dm 2 


h' 2 h 2 


d cos 6 d(p 





(1) 



where hi{h'^) and H(H') are indices labeling the helicity states of quark and nucleon, and 
Hi{H[) labeling the helicity state of resonance the (a, p). See Fig. |l[ In order to incorporate 
the final state interaction, we have separated the q — > 7r + 7r~ fragmentation process into two 
steps. First, the quark fragments into the resonance (a, p) , then the resonance decays into 
two pions, as shown at the top part of the Fig. [I]. 

We first discuss two-meson fragmentation, first examined in Ref. @. Here we introduce 
only those pieces necessary to describe s-p interference in 7r + 7r~ production. A full account 
of these fragmentation functions will be given in Ref. [11]. A two- meson fragmentation 
function can be defined by a natural generalization of the single particle case. Using the 
light-cone formalism of Collins and Soper |13|], the following replacement suffices, 



\hX) out out (hX\ 



\n + H X)out out(7T + 7T X\. 



(2) 



The resulting two meson fragmentation function depends on the momentum fraction of each 
pion, zi, Z2, the irn invariant mass, m, and the angles G and <fi. The subscript "out" places 
outgoing wave boundary conditions on the tcttX state. Two types of final state interactions 
can generate a non-trivial phase: i) those between the two pions, and ii) those between the 
pions and the hadronic state X. We ignore the latter because we expect the phase to average 
to zero when the sum on X is performed - - |/T + 7r~X) out — > |(7r + 7r~) out X). Furthermore, if 
the two-pion system is well approximated by a single resonance, then the resonance phase 
cancels in the product |(7r + 7r~) out X) ((7r + 7r~) out X|. This leaves only the interference between 
two partial waves as a potential source of an asymmetry. The final state phase of the two- 



pion system is determined by the rnr T- matrix ||14|| . We separate out the phase for later 
consideration and analyze the (real) p-a interference fragmentation function as if the two 
particles were stable. 

The s-p interference fragmentation function describes the emission of a p(o~) with helicity 
Hi from a quark of helicity h 2 , followed by absorption of er(p), with helicity H[ forming a 
quark of helicity h' 2 . Conservation of angular momentum along the £3 axis requires 



H x + h' 2 = H[ + h 2 . 
Parity and time reversal restrict the number of independent components of Xi: 



M 



sp(ps) 
H[Hi,h 2 h' 2 



M 



sp(ps) 

~H' 1 -H 1 -h 2 -h , 2 



(parity) 



Jvl H[H 1 ,h 2 h' 2 



M P H lHi>h > 2h2 (T- reversal) 



(3) 

(4) 
(5) 
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Note that Eq. ([5]) holds only after the T-reversal violating final state interaction between two 
pions is separated out. After these symmetry restrictions, only two independent components 



MZ, ++ = A< ++ = MZ,- = A<__ oc q, , (6) 
MZ, + - = MZ,- + = Mti,- + = M p _ s 10 , + _ ex 6q T , (7) 

and they can be identified with two novel interference fragmentation functions, q n Sq n where 
the subscript / stands for interference. Here, to preserve clarity, the flavor, Q 2 , and z have 
been suppressed. The helicity ±~ states of quarks are denoted ±, respectively. Hermiticity 
and time reversal invariance guarantee q I and 8q 1 are real. 

From Eq. (0) it is clear that the interference fragmentation function, <5g 7 , is associated 
with quark helicity flip and is therefore chiral-odd. It is this feature that enables us to access 
the chiral-odd quark transversity distribution in DIS. 

Encoding this information into a double density matrix notation, we define 

d 2 M 

Ao(m 2 ){I®fioq 1 (z) + (a + ®fj- + a-®fj + )5q I (z)}A* 1 (m 2 ) 



dz dm 2 

+A 1 (m 2 ) {I ® 770 q z {z) + (<r_ ® V+ + °+ ® V-) A* Q (m 2 ) , (8) 

where a± = (<7i ± io2)j2 with {a{\ the usual Pauli matrices. The ??'s are 4x4 matrices 
in (a, p) helicity space with nonzero elements only in the first column, and the 77's are the 
transpose matrices (7^0 = Vo:V+ = V T iV~ = V+)y with the first rows (0,0,1,0), (0,0,0,1), 
and (0, 1,0,0) for fj , fj + , and f/_, respectively. The explicit definition of the fragmentation 
functions will be given in Ref. [|11[ . 



The final state interactions between the two pions are included explicitly in 

A (m 2 ) = -ismS e i5 ° , Ax(m 2 ) = -isin^ 1 , (9) 

where 5q and 8\ are the strong interaction irn phase shifts. Here we have suppressed the m 2 
dependence of the phase shifts for simplicity. The decay process, (a, p) — > irir, can be easily 
calculated and encoded into the helicity matrix formalism. The result for the interference 
part is 



d 2 V V6 



d cos Odd) 87r 2 m 



sin G 



te 



i<t> ( ?7 _ _ fj_) + ie i<fi ( v+ _ n+ ) - V2 cot 6 (rlo + 77o)l • (10) 



Here we have adopted the customary conventions for the p polarization vectors, e± = =F(ei ± 
ie 2 )/\/2 and e* — &3 m its rest frame with e's the unit vectors. 

In the double density matrix notation, the quark distribution function T can be expressed 
as H 

T = ]-q{x) ^Aq(x) a 3 ® a 3 + ^Sq{x) (cr+ ® cr_ + a_ ® a+) , (11) 

where the first matrix in the direct product is in the nucleon helicity space and the second 
in the quark helicity space. Here q, Aq, and 5q are the spin average, helicity difference, 
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and transversity distribution functions, respectively, and their dependence on Q 2 has been 
suppressed. 

The hard partonic process of interest here is essentially the forward virtual Compton 
scattering as shown in the middle of Fig. [|. For an unpolarized electron beam, the resulting 
cross section is M 



d 2 a(eq — > e'q') e*e; 



p 4 ^ 2 ri l ft _ 



dx dy 8nQ 2 



l + (l-y) 2 lT ^ T , _ . , 2(1 -y) 



2y 



7 <g) J + cr 3 (g cr 3 ) H (cr + ® ct_ + ct_ (g) 



(12) 



in the collinear approximation. Here e q is the charge fraction carried by a quark. We have 
integrated out the azimuthal angle of the scattering plane. 

Combining all the ingredients together, and integrating over to eliminate the q T de- 
pendence, we obtain the transversity dependent part of the cross section for the production 
of two pions (kaons) in the current fragmentation region for unpolarized electrons incident 
on a transversely polarized nucleon as follows 

d 5 a± 7T e 4 1 - y /- 

VD COS0 



dx dy dz dm 2 dcj) 2 327r 3 Q 2 m y 

x sin5 sin5i sin (S - Si) ^e 2 a Sq a (x) 6q*(z) . (13) 

a 

Here the sum over a covers all quark and antiquark flavors. 

An asymmetry is obtained by dividing out the polarization independent cross section, 

da ± -da T 7T VQ(l~y) , . x . x . , x X s 

A± T = -j — i — = -- A -r—r, vo cos0 sm S sin dx sm (S -Si) 

da± + dcj 4 1 + (1 — y) 1 

E a e 2 a Sq a (x)Sq%z) 



Ea elq a (x) [sin 2 5 qg(z) + sin 2 Siqf(z) 

where go and qi are spin-average fragmentation functions for the a and p resonances, re- 
spectively. This asymmetry can be measured either by flipping the target transverse spin 
or by binning events according to the sign of the crucial azimuthal angle <fi. Note that 
this asymmetry only requires a transversely polarized nucleon target, but not a polarized 
electron beam. 

The flavor content of the asymmetry A±_t can be revealed by using isospin symmetry and 
charge conjugation restrictions. For 7r + 7r~ production, isospin symmetry gives Sii I = —Sdj 
and Ssj = 0. Charge conjugation implies Sq a = —Sq . Thus there is only one independent 
interference fragmentation function for 7r + 7r~ production, and it may be factored out of the 
asymmetry, J2 a e !Sq a Sq® = [4/9 (<kt — Su) — l/9(Sd — Sd)]Su r Similar application of isospin 
symmetry and charge conjugation to the p and o fragmentation functions that appear in 
the denominator of Eq. (|14|) leads to a reduction in the number of independent functions: 

Ui = di = Ui = di and = Sj for i = {0, 1}. For other systems the situation is more 
complicated due to the relaxation of Bose symmetry restrictions. For example, for KK 
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system, 5qf = —dq 1 still holds, but 5uj, Sdj, and Ss n are in general independent. We also 
note that application of the Schwartz inequality puts an upper bound on the interference 
fragmentation function, Sq^ < 4g <?i/3 for each flavor. 

Finally, a few comments can be made about our results. First, the final state phase 
generated by the s-p interference is crucial to this analysis. If the data are not kept differ- 
ential in enough kinematic variables, the effect will almost certainly average to zero. We are 
particularly concerned about the two-meson invariant mass, m, where we can see explicitly 
that the interference averages to zero over the p as shown in Fig. |2|. Second, the transver- 
sity distribution is multiplied by the fragmentation function Sq T . Note that the transversity 
distribution always appears in a product of two soft QCD functions due to its chiral-odd 
nature. In order to disentangle the transversity distribution from the asymmetry, one may 
invoke the process e + e" — > (7r + 7r~X)(7r + 7r~X) to measure 8q n or use QCD inspired models 



to estimate it 11 



To summarize, we have introduced twist-two interference quark fragmentation functions 
in helicity density matrix formalism and shown how the nucleon's transversity distribution 
can be probed through the final state interaction between two mesons (tt + tt~, ttK or KK) 
produced in the current fragmentation region in deep inelastic scattering on a transversely 
polarized nucleon. The technique developed in this Letter can also be applied to other 
processes. Straightforward applications include the longitudinally polarized nucleon, and 
e + e _ — > (7r + 7r~X)(7r + 7r~X). A somewhat more complicated extension can be made to 
two-meson production in single polarized nucleon-nucleon collisions — pp± — > tt + tt~X, etc. 
These applications will be presented in future publications. 



We would like to thank John Collins and Xiangdong Ji for helpful conversations relating 
to this subject. 
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FIG. 1. Hard scattering diagram for it + it~(KK) production in the current fragmentation 
region of electron scattering from a target nucleon. In perturbative QCD the diagram (from 
bottom to top) factors into the products of distribution function, hard scattering, fragmentation 
function, and final state interaction. Helicity density matrix labels are shown explicitly. 



to 




o 




52. 


0.6 


E 




m 




vo 

c 


0.4 


w 




o 

to 




C 


0.2 


CO 











0.2 




0.4 







.!,,,!,,,!,, 


, I 7 ,,!,,, I 



0.5 0.6 0.7 0.8 0.9 1 

m(GeV) 



FIG. 2. The factor, sin<5o sin<5i sin((5o — 6\), as a function of the invariant mass m of 
two-pion system. The data on -kit phase shifts are taken from Ref. ITU 
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